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TECHNICAL NOTE D-92} 

AN AUTOMATIC TERMINAL GUIDANCE SYSTEM FOR RENDEZVOUS 
WITH A SATELLITE 
By Terrance M. Carney 

SUMMARY 


This study includes a consideration of the design philosophy for 
an automatic terminal guidance system, a derivation of guidance equa- 
tions required, and an outline of the general type of instrumentation 
necessary to provide the essential information. A control system for a 
sample vehicle is analyzed. 

A representative case, rendezvous with a satellite in circular 
orbit at 400 nautical miles, was examined. Terminal- stage nominal 
burning times of 200 and 400 seconds were used. 

For the 200-second case, initial errors in circumferential dis- 
placement of ±25,000 feet, in radial displacement of 7 ,000 to -9,000 feet, 
and in lateral displacement of ±20,000 feet were within the capabilities 
of the system. Velocity errors of 300 to -400 ft/sec in the circumfer- 
ential direction, 180 to -200 ft/sec in the radial direction, and veloc- 
ity offsets of at least 2° (±800 ft/sec) in the lateral direction could 
also be handled. The 400- second case was capable of correcting larger 
errors, but limits were not determined. 

The dependence of required characteristic velocity on initial 
errors was determined and it was found that increases over the nominal 
terminal- stage characteristic velocity of the order of 15 percent 
covered most of the previously mentioned in-plane errors. The require- 
ments were more severe for cases with lateral velocity offsets. A 
simplified set of guidance equations was tested and produced only slight 
variations in performance. 

Overall velocity requirements and mass ratios were determined for 
terminal- stage burning times of 100, 200, and 400 seconds and for 

a range of transfer angles by using exact calculations for the terminal 
stage and an impulsive launching velocity. These results indicated that 
the shortest burning time consistent with the launch guidance errors 
expected gave the best mass ratio. 
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INTRODUCTION 


One of the many problems engendered by space operations is that of 
physical comminl cation between orbiting space stations and the earth. 

The most straightforward means of accomplishing this ccranu ni cation is 
by direct launch of a vehicle from earth and subsequent guidance into 
coincidence in velocity and position with the orbiting station. This 
problem and the corollary problems of bringing orbiting vehicles together 
in space (orbital transfer) and of matching position but not the veloc- 
ity of the station (hard rendezvous) have received much attention in the 
recent literature. References 1 to 6 are representative of work being 
done in this area. 

The overall problem of launch and guidance of a vehicle to rendez- 
vous is one of large scale and is best handled by tearing it into several 
phases. A reasonable partitioning is into launch, midcourse guidance, 
terminal guidance, and docking. There is a strong interaction between 
these phases which must always be checked, but for first considerations 
this separation is convenient. 

This paper will be concerned with the terminal guidance stage, which 
is here defined to conmence with onboard sensor acquisition of the target 
vehicle prior to firing of the final stage and to extend to the beginning 
of the docking phase. Both automatic and piloted systems have previously 
been proposed for this task (refB. 1, 3> and 4); only the automatic 
approach will be treated here. 

Included in this study is a consideration of the design philosophy 
for a term inal guidance system, a derivation of guidance equations 
required, and suggestions for the general type of instrumentation 
required to provide the necessary information. 

The following steps were taken in the design and examination of 
this system: 

1. Closed- form solutions were written for simplified equations of 
motion. 


2. These solutions were compared with exact numerical sdlutions 
to determine that no gross errors had been introduced by the 
simplifications . 

3. Guidance relations were formed from the analytical solutions 
by substitution of measured variables. 

4. A control technique for steering the rendezvous vehicle toward 
the trajectory defined by the guidance relations was devised and 
analyzed. 
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5. A typical case was examined by programing simplified equations 
of motion for the IBM type 704 electronic data processing machine and 
testing the proposed system for a variety of initial errors with respect 
to a nominal initial aim point for terminal-stage burning times of 200 
and 4-00 seconds . 

6. The penalty associated with using extended terminal- stage 
burning times was determined by calculating the overall mass ratio for 
various terminal- stage burning times plus an associated impulsive launch. 


SYMBOLS 


The axis system and conventions used in this report are illustrated 
in figures 1, 2, and 3 . The symbols used are defined as follows: 

A constant in guidance equations 

c effective rocket exhaust velocity, ft/sec 

d x displacement of attitude control thruster from center of 

gravity along x body axis, ft 

e unit vector 

f thrust per unit mass of commuter vehicle, lb/slug 

thrust per unit mass along line-of-sight projection in 
XY-plane , lb/ slug 

thrust per unit mass normal to line of sight in XY-plane, 
lb/ slug 

acceleration due to gravity, ft/sec^ 
apparent gravitational acceleration, ft/sec^ 
unit vector along line-of-sight projection in XY-plane 
moments of inertia about x, y, z body axes, respectively, 

slug-ft^ 

specific impulse, sec 

unit vector normal to line-of-sight projection in XY-plane 
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My,M z 


H 

R 

s 

t 

t g 

At 

T 

u,v,v 

V 

AV 


x,y,z 


X,Y,Z 


unit vector parallel to Z-axis 

gain constants in longitudinal controller 

navigation constant in lateral controller , sec 

attitude control thrust slope, lb/radian 
unit vector along V c ^ s 
mass, slugs 

control moments about y and z body axes, respectively, 
ft -lb 

angular rates about x, y, and z body axes, respectively, 
radians/ sec 

performance function 

radial distance from center of earth, ft 
control sensitivity 
time, sec 

time remaining to rendezvous from reference time , sec 

burning time, sec 
thrust, lb 

velocities along x, y, z body axes, respectively, ft/sec 
total velocity, ft/ sec 
characteristic velocity, ft/sec 

body axes with origin at commuter center of gravity and 
displacements along these axes 

rectangular axes with origin at space station at t = 0 
and displacements along these axes 


a 


elevation angle, radians 
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p 

angle of lateral offset in XY-plane, radians 

7 

gravitational constant , ft^/sec^ 

8 

attitude control valve displacement, radians 

A 

error 

e 

error signal 

£ 

ratio of damping to critical damping 

,T ^z 

total impulse consumed by axial, yaw control, and pitch 
control thrusters, lb-sec 

e,i,0 

Euler angles, radians 

A 

Laplace operator 

n 

mass ratio, ratio of mass at launch to final mass 

p 

line-of-sight range, ft 

cr 

lateral velocity offset angle, deg 

T 

characteristic time, sec 

* 

yaw angle with respect to line of sight, radians 

0) 

natural frequency, radians/sec 

n 

angle of transfer, deg 

Qjp 

angle of transfer from launch to initiation of terminal 
stage, deg 

Subscripts : 


c 

commuter 

cZ 

closed loop 

com 

command 

cs 

control system 

e 

earth 
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f at rendezvous 

i at initiation of terminal stage 

I inertial reference 

Is line of sight 

L yaw servo gain 

meas measured 

n normal 

nom nominal initial condition 

0 at reference t = 0 

01 open loop 

p pitch servo gain 

q pitch rate gain 

r yaw rate gain 

R radial 

s space station 

t circumferential 

x,y, z along or about x-, y-, or z-axes, respectively 

X, Y, Z along or about X-, Y-, or Z-axes, respectively 

0 pitch angle gain 

t yaw angle gain 

A dot over a variable indicates a derivative with respect to time; 

two dots indicate a second derivative with respect to time. 

A bar over a variable indicates a vector. 
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ANALYSIS 


System Philosophy 

The ideal terminal guidance system is one that follows a nominal 
path which is optimum from the fuel- consumption standpoint, allows 
correction of gross errors in initial conditions without instability 
or excessive additional fuel use, and employs highly reliable instrumen- 
tation in minimum quantity. It is clear that these criteria must be 
compromised to produce a system which will best satisfy them in combina- 
tion. This can be illustrated by the contrast between the optimum-fuel- 
consumption path, which consists of an impulsive boost to orbital veloc- 
ity at apogee of the launch trajectory, and the best path for error 
correction, which consists of a very long low-acceleration burning 
period during which the thrust vector can be controlled. The impulsive 
boost admits of no displacement-error correction capability, while the 
low-acceleration path is rather inefficient from the standpoint of fuel 
consumption. Again, a complete inertial reference offers the possibility 
of very sophisticated steering techniques, but introduces drift error 
problems as well as weight and volume penalties, while simple propor- 
tional navigation requires no fixed reference but steers a path which 
is not the most economical of fuel. 

The ter minal guidance system of this paper was developed around a 
variable- thrust rocket motor, since control of both direction and magni- 
tude of the t hrus t vector provides maximum flexibility in correction 
of errors. r n» a philosophy adopted in the design of the system was that 
a path requiring the least steering of the cammiter vehicle should be 
used, that the thrust demanded during flight should be a smooth, slowly 
varying function, and that the simplest reference possible should be 
employed. 


Mass Particle Solution 

In order to meet the requirements of minimum steering and a smooth 
thrust function, the nominal flight path was defined as one where con- 
stant thrust was maintained in a constant direction with respect to a 
local earth reference. The simplest case, rendezvous with a station 
in circ ular orbit, was chosen because it is most amenable to closed- 
form solutions. 

A solution to simplified equations of motion of a particle with 
constant thrust and linear mass variation with time was undertaken to 
study this f light , path. These equations were written in a rectangular 
axis system (fig. l) where it was assumed that the motion of the target 
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or space station is linear along the X-axis. It was further assumed, 
that the force of gravity is constant in the neighborhood of the space 
station, and centrifugal effects were introduced by modifying the 
gravity term. The equations are: 


Z 


X = 


Y = 


mo 


+ mt R 


T X 

(1) 

i Hq + mt 


Ty 

(2) 

iDq + mt 


* - *0 

(3) 


In reference 7 it is shown that thrusting along the velocity vector of 
the vehicle is an efficient method of gaining velocity. In the terminal 
stage of rendezvous with stations in nearly circular orbit the velocity 
vector lies substantially in the local horizontal plane. Since the 
local horizontal is a convenient reference and some complication is 
involved in determining the true velocity direction, thrusting in the 
horizontal plane provides a good approximation to the efficient course 
and is relatively simple to mechanize. Therefore Ty and T^ were 

set equal to zero and the relationship between Ty and m written 

T x = -me (4) 

where c is the effective rocket exhaust velocity. A solution in the 
XZ- plane is obtained first. The required end conditions are 

X(t f ) = V s 
X(tf) - v s t f = 0 
Z(tf) = Z(t f ) = 0 


where t f is defined as the time at rendezvous. Equations (l) and (3) 
were integrated and these required end conditions applied to obtain the 
following closed- form solutions for the required geometrical initial 
conditions in terms of burning rate m and remaining flight time t g . 

io = v s + C log e (l + ± t g ) (5) 
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Inspection of these four equations reveals that fixing any two of 
the six parameters Xq, X q , £ 0 , Z 0 , t g , 8X1(1 m defines a unique 

solution. The time of flight can then be selected and the initial 
circumferential velocity Xo can be determined from launching trajec- 
tory considerations to calculate nominal initial conditions. In the 
real case, the four geometrical variables will be fixed at the time of 
radar acquisition. The procedure proposed is that Xq and X Q be 
considered the prime variables and. that measurements of these be used 
to compute continuously the Z Q , Z Q , and thrust magnitude T x required. 

Appropriate directional control of the thrust vector will then be exer- 
cised to drive the measured Z and Z to the required values. When 
this condition is satisfied, the commuter vehicle will be on a flight 
path requiring no further steering and using constant thrust in the 
local horizontal plane to rendezvous with the target station. The rela- 
tions required, then, are Z Q , Z Q , and T x in terms of Xq and o' 

Solving equations (40 and ( 5 ) for T^m gives 
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Express t g by using equations (5) and (6) and substituting a series 

representation for the exponential, and manipulate equations (7 ) » (®)> 
and (9) to yield the following expressions: 




Harms above order [(*• - *o)A] 2 have been dropped since this cor- 
responds to order 10 - ^ for the rendezvous situations considered in this 
paper. 

In order to e valuat e the assumptions used to obtain the prece din g 
solutions, exact solutions for coplanar rendezvous with a station in a 
HOO-naut ical-mile , circular orbit with 200- and HOO- second terminal- 
stage burning t ime s were calculated numerically . The exact solutions 
were for c ons tant thrust in the horizontal plane with the constraint 
that the initial circumferential velocity should match the analytical 
solution. The equations and techniques used in arriving at these solu- 
tions are outlined in appendix A. Die geometric errors introduced by 
the assumptions are illustrated in figure 4. For the 200- second case 
the mmfiMiM - radial velocity 4 difference is 8 ft/sec and the maximum 
radial displacement Z difference is 70 feet, notice that the cor- 
respondence Improv es near rendezvous. Die agreement for the 400- second 
case is not as good, showing a maximum radial velocity difference of 
98 ft /sec and maximum radial displacement of 1,000 feet. In operation 
the system continually computes the analytic solution, and as the time 
to rendezvous decreases the analytic solution will approach the cor- 
responding exac t solution until near rendezvous they will he effectively 
coincident. 

Hip errors in thrust required were also studied. Satisfaction of 
the velocity constraint mentioned previously required thrusts respec- 
tively 2-9 *»nd 12.1 percent higher for the exact 200- and 400- se cond- 
huming-time cases than for the correspond 1 ng analytic solutions. Dils 
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thrust error Is attributed to the neglect of coupling of the radial 
motion in the equation describing circumferential velocity gain in the 
simplified system. The closed- loop guidance would tend to nullify 
these errors since it would predict successively more accurate thrusts 
as the vehicles closed. 

The effect of this thrust discrepancy on computed fuel consumption 
was also determined. Figure 5 is a plot of the ratio of characteristic 
velocity for the exact case to the initial circumferential velocity of 
the commuter relative to the station against the initial circumferential 
velocity difference for terminal- stage burning times of 100, 200, 300, 
and 400 seconds. The expression used for calculating characteristic 
velocity is: 

™ - sAp l°Se(^) 

' c l°8e(^) <«) 

Inspection of equation (5) shows that the characteristic velocity cal- 
culated in the rectangularized solution will be the initial circumfer- 
ential velocity relative to the station. The regularity of the varia- 
tion of the ratio plotted in figure 5 indicates that simple corrections 
can be applied to fuel consumptions calculated in the rectangularized 
system to make them significant. 

The next consideration was guidance in the XY- plane. The fact 
that the Y-dynamics were not coupled to the XZ-system led to the con- 
clusion that offsets in Y and Y could be handled independently of 
the in-plane guidance, .parti cularly in the case where a relatively long 
burning t ime is used and rapid response is not demanded of the vehicle. 
Since maneuvers in the XY-plane are not coupled to gravity (to first 
order) the path which is followed does not strongly affect the effi- 
ciency of the system as long as the path is smooth and does not require 
large or frequent accelerations transverse to the X-axis. Wrigley, in 
reference 5, has described an application of the proportional navigation 
technique to satellite rendezvous. The advantage of such a system is 
that it does not require a measurement of bearing angle which depends 
on an inertial reference. Only information about the range and rate of 
change of range to the target and the rate of change of bearing of the 
line of sight is needed. This technique, when restricted to the 
XY-plane, was believed to be satisfactory for control of Y- and 
Y-offsets. Development of this system for the current application is 
carried out in appendixes B and C. The control equation yielded is: 
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Pp 


x 

m 


(14) 


where ? is the yaw tilt angle of the thrust vector from the line of 
sight in the horizontal plane. 


It is clear that the elliptic case is more difficult to. study 
analytically than is rendezvous with a station in circular orbit. It 
is submitted, however, that for small eccentricities this problem can 
be handled by biasing the XY reference plane to a direction parallel 
to the inclination of the desired velocity vector at rendezvous to the 
local horizontal. This inclination can be determined as a function of 
the known orbital elements of the station and the flight time plus the 
predicted remaining flight time t g of the commuter vehicle. 


Instrumentation 

In order to obtain the accuracy required in the terminal stage of 
rendezvous it is best to measure the guidance information in relative 
coordinates from onboard the commuter or station. Measurement of the 
required quantities from the earth would introduce a command delay time, 
as well as inaccuracies from long ranges and extra coordinate 
transformations . 

The particle solutions defined the quantities which must be meas- 
ured for guidance purposes. A general plan for providing these meas- 
urements was evolved and is illustrated in figure 2. The local vertical 
reference is provided by a horizon scanner. Radar or a similar technique 
is used to measure range and range rate. The elevation of the line of 
sight from the local horizontal must be measured. The rates of slewing 
of the line of sight with respect to inertial space in the XY-plane 
and with respect to the local horizontal in elevation are required. In 
the cases studied, rates of the order of a milliradian/second were 
encountered. To provide sufficient precision to control these rates, 
measurement accuracies of the order of 0.1 milliradian/second or better 
are required. If conventional means of rate measurement, such as 
mounting gyros on a radar antenna or differentiating successive angular 
measurements, will not yield this accuracy, an auxiliary system, such 
as an optical tracker, may be necessary. 

Difficulties arising from the noise properties and pulsed charac- 
ter of the radar information and other instrument errors are not con- 
sidered in this paper. It is believed that these effects can be 
counteracted without changing the character of the system, and that 
most of the measurement errors will be in part compensated by the 
closed- loop configuration. 
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The particular mission plan will determine which vehicle carries 
the sensing equipment. An instrumented commuter would he required for 
rendezvous with passive stations, while rendezvous of several commuters 
with one space station might make location of the sensors at the station 
more economical. In either case, the guidance technique described can 
be used. 


L 

1 

5 

2 
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Guidance Equations 

Once the parameters to be controlled have been expressed and the 
measured variables have been defined, it is necessary to develop rela- 
tions for these parameters in terms of the actual sensor inputs. It is 
also desirable to find relations which are relatively simple to compute 
with onboard equipment. In order to do this, measured (line of sight) 
variables are introduced into the mass particle solutions for T x , Z, 

Z, and the lateral command quantity ty. The trigonometric functions, 
involved have been replaced by their small-angle approximations. This 
procedure is carried out in appendix C. The resulting expressions are: 


^com ” 




(15) 


where 



(16) 
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m /com 
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(17) 
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Y com 0 t x 

HT 


(18) 


where equation (l8) is equation (l 4 ) repeated. 


Note that the expression for axial acceleration is similar to that 
used by Spradlin in reference 1 . 


Although these relations could be generated by a suitable analog 
computer, the large range covered by the variables would seem to make 
digital computation more appropriate. The sampled data form of the 
radar output also lends itself to this computation technique, since a 
digital computer could be used for data smoothing. 


Commuter Stability Analysis 

Determination of the guidance commands led to the next requirement, 
a steering system to compare the commands to measured values and to 
correct the observed errors. The control system proposed operates by 
changing the commuter vehicle's attitude to provide components of the 
axial thrust in the proper direction to correct error signals generated 
by the guidance computer. The magnitude of the main engine thrust is 
determined by equation (17) without compensation for pitch or yaw angles. 
Attitude control is provided by a set of four thrusters at the nose of 
the vehicle. (See fig. 3.) Roll stabilization of the vehicle with 
respect to the local horizontal is assumed. This mode of control was 
chosen to avoid the complication of gimballing the main rocket, and 
because it was slightly less complex than a pure torque attitude con- 
troller, which requires four additional thrusters. 

An elementary analysis of the dynamics of the controlled vehicle 
was carried out by means of the root- locus technique to determine means 
of stabilizing the vehicle and system gains which provide satisfactory 
response. Block diagrams of the longitudinal and lateral controllers 
investigated appear in figure 6 and the analysis appears in appendix D. 

It was necessary to compromise between requiring such low system gains 
that no initial errors could upset the vehicle and providing high enough 
gains to insure rapid tracking near the end of the terminal stage. Nor- 
mally, this is best achieved by use of variable gains or nonlinear com- 
pensation. The generalized system considered- in this paper used linear 
gains with limits on the error signals. Hiese limits were arbitrarily 
set so that no tilt in pitch greater than 0.6 radian or tilt in yaw 
greater than 0. 5 radian could be called for. The maximum force-to-mass 
ratio of the attitude-control thrusters was also limited to 0.25 pound 
per slug to maintain a reasonable sizing relative to the main motor. 
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These limits were checked only to determine their effect on stability, 
and no attempt was made to attain the optional limited system. 


RESULTS AMD DISCUSSION 


In order to demonstrate the flexibility of this terminal guidance 
concept and to detect any inherent problems in its usage, a particular 
rendezvous case was tested. The case chosen was terminal rendezvous 
with a station in 400- nautical-mile circular orbit. Nominal initial 
conditions were based on an 82.5° transfer from impulsive launch. Ref- 
erence 6 points out that a transfer of this order offers advantages in 
frequency of possible rendezvous, and the orbital altitude was selected 
as representative for a space station. A commuter vehicle with a mass 
of J00 slugs was assumed with reasonable associated physical 
characteristics. 

The equations of motion in body axes were programed for numerical 
solution on the IBM type 704 electronic data processing machine, and 
solutions were carried out in terms of the recta n gu l a r reference axes 
previously discussed, lhe equations and physical characteristics of 
the system studied are included in appendix E. This scheme was employed 
rather than exact equations to simplify the programing and because the 
previously noted investigation had shown the geometric errors introduced 
to be reasonable and had indicated that a simple linear correction (see 
fig* 5) applied to the fuel consumption measured in the simplified 
system would make it correspond well to exact values. Moreover, most 
of the data computed were evaluated on a comparative basis so that 
deviations of the nominal values from the exact cases did not affect 
any conclusions. 

Nominal terminal- stage burning times of 200 and kOO seconds were 
studied, with emphasis on the 200-second case. Although introduction 
of errors change s the burning time, all cases related to a given 
nnwrinai burning time will be referred to by the nominal time for con- 
venience. For the 200- second case, the acceptable XZ-plane error 
perturbations in velocity and displacement from the nominal initial 
conditions were determined. Two lateral velocity offsets, ^corresponding 
to the conditions where the commuter velocity was 1° and 2 out of the 
orbital plane of the station, and various side displacements were also 
studied. The relation of initial side displacement to initial lateral 
velocity offset for minimum fuel consumption was determined. Comparable 
coplanar and offset cases were run for the ^JOO- second burning time. 
Nominal initial conditions for these cases appear in table I. The fuel 
consumption was expressed as a characteristic terminal velocity require- 
ment for 1 these causes to determine variations with Initial errors. 

A simplified version of the guidance scheme was tested for the 200-second 
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case "to determine efficiency loss due to leaving out some of the terms. 
Finally, the overall mass ratio variation with terminal burning time 
was determined by including an impulsive launch. 


In-Plane Capability 

The initial error correction capability of the system was checked 
by the simple technique of perturbing each initial condition with suc- 
cessively larger errors until the system was no longer able to correct 
itself. This is not representative of the real case, where errors 
vill generally occur in combination, but does indicate limits without 
the difficulty of determining the most probable error combinations. 

Figures 7 and 8 illustrate trajectories and thrust profiles for 
nominal initial conditions and limiting errors for the 200-second 
burning time. Figure 7 shows the nominal trajectory and four error 
perturbations in initial displacement, while figure 8 shows the nominal 
trajectory and four error perturbations in initial velocity. Veloci- 
ties and displacements shown are measured relative to the space station. 
For this case, the range of acceptable circumferential errors is 
25,000 to -25,000 feet in displacement and 500 to -400 ft/sec in veloc- 
ity. The radial-error range is somewhat more restrictive, being 7,000 
to -9,000 feet in displacement and 180 to -200 ft/sec in velocity. The 
400-second case was checked for the same extremes as the 200-second 
case. This case is theoretically capable of correcting larger velocity 
and displacement errors during its longer operating time, but its limits 
were not determined. 

Inspection of the thrust time histories for the various errors in 
figures 7(a) and 8(a) reveals that the ratio of maximum to minimum 
thrust required is 3«25* This is well within the 5-1 range of modula- 
tion reported in reference 8 for an operating variable -thrust rocket 
motor. The smooth character of the thrust profiles indicates that the 
dynamic response demanded of the throttling system is relatively slow. 

Figures 7(b) and 8(b) are time histories of the last 10 seconds 
before rendezvous. These are included to show the magnitude of the 
closing velocities and displacements in more detail, and to indicate 
the size of residual errors to be expected. No residual velocities 
greater than 1 ft/sec or displacements greater than 5 feet were encoun- 
tered. However, the large closing velocities shortly before rendezvous 
indicate that a bias should be introduced to avoid collision of the 
vehicles due to measurement errors. If this bias were used, the com- 
muter would rendezvous with a point a short distance ahead of the sta- 
tion, and the docking system could then perform vernier corrections and 
direct coupling of the vehicles. 
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The limiting factor in these cases proved to be control-system 
stability. For large inputs, the limited system behaved as an on-off 
controller and, with the control sensitivities used, the system tended 
to switch too frequently and did not succeed in nullifying the Z and 
2, errors before the space station closed on the commuter for initial 
errors outside the stated bounds. In reference 9 j Schmidt outlines a 
technique for designing a nonlinear compensator which operates on the 
error signal to force the limited system to follow the optimum switching 
solution. For a particular system, application of this procedure should 
considerably improve the response of the controller and extend its 
error-handling capability. 

Attitude -control- system impulse requirements were calculated for 
all cases and never exceed 1,000 pound-seconds for the proportional 
system used. Variations at this level are insignificant compared to 
variations in fuel consumption of the main motor. 


Cross-Plane Capability 

The ability of the system to handle velocity offsets of 1° and 2° 
was tested for a range of initial displacements. Figure 9 shows the 
1° offset case for a nominal burning time of 200 seconds. Values of 
Kp, the navigation constant, from 2.0 to 7-0 were tested. Marginal 

stability was exhibited for K£ = 2.0. Curves are plotted for Kp = 3-0 

and 7.0 to illustrate variations due to gain change. Comparison of 
fuel consumption for stable values indicated that = 3.0 was most 

economical, and no significant difference between values could be seen 
from a stability standpoint. For any cases tested with cross-plane 
errors, the initial conditions were adjusted so that the longitudinal 
guidance system detected no initial errors. HLere was little apparent 
coupling of the motion in the two planes during any of the cases for 
errors within the system limitations. 


It was determined from inspection of fuel- consumption variation 
that the opt imum initial condition for the offset case was that the 
initial angular rate of the line of sight be zero, that is, 


Po ~ 


Vo * *0 ( V B - Xp) 

X 0 2 ♦ Y„ 2 


- = 0 


(19) 


or 



Xq - V s 


( 20 ) 
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Inspection of the thrust time histories in figure 9 shows that the 
optimum initial condition (Yj_ = 39>537) demands nearly constant thrust 
throughout. Comparison with the nominal case for coplanar rendezvous 
in figures 7 and 8 shows that a nearly constant increment of additional 
thrust amounting to about 10 percent is required to correct the velocity 
offset. 

If the optimum initial condition is defined as nominal, the system 
is capable of handling initial Y-displacements of ±20,000 feet for 
initial velocity offsets of 0°, 1° (Yo = -408.2 ft/sec), or 2° 

(Yo = -8l6.5 ft/sec). Time histories of the final 10 seconds before 
rendezvous appear in figure 10. These are for the 2° lateral velocity 
offset case, which yields larger residual errors than the 1° offset 
case. Final velocities as large as 9 ft/sec and displacements up to 
8 feet occurred for the most severe errors. This result further indi- 
cates that a bias should be introduced to prevent collision of the 
vehicles. The error correction capability in the cross-plane case is 
determined by the maximum available lateral thrust of the system at the 
prescribed- limit yaw tilt of 0.5 radian. Methods of increasing this 
tolerance were not sought, since it was already rather broad. Larger 
velocity offsets were not tested, and again the ^-00-second-burning-time 
case should yield larger initial tolerances but was not tested. 


Characteristic Velocity Variation With Initial Errors 

In order to illustrate the dependence of fuel consumption on 
initial errors, the characteristic velocities required for the condi- 
tions tested are plotted in figures ll(a) and ll(b). The characteristic 
velocities were calculated by using equation ( 13 ) and were based on the 
mass ratios from the computer runs. No corrections were introduced. 

This form is intended to generalize the results as much as possible 
since it removes direct dependence on specific impulse. Although the 
particular effective exhaust velocity used (c = 10,000) will affect 
the results due to coupling of the rate of change of mass in the gravity 
field, mass ratios based on the included data and scaled to specific 
impulses in this neighborhood should be good approximations. 

The variation of required characteristic velocity with errors 
plotted in figure ll(a) shows that the in-plane errors generate almost 
linear slopes, with the exception of the variation with initial vertical 
velocity which is roughly parabolic. If it can be guaranteed that the 
circumferential velocity will be low by no more than 150 ft/sec, the 
maximum errors require 15 percent additional characteristic velocity 
for the 200-second-burning-time case. For the 400-second case, the 
same errors require 13.2 percent additional characteristic velocity, 
and since this case is based on a lower nominal value, there is an 
additional savings in magnitude of the added velocity required. 
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For the cases with lateral velocity offsets the sharp increase in 
required characteristic velocity which occurs with increasing offset 
angle is illustrated in figure ll(b) . The variation of required velocity 
with displacement errors is also large, particularly for the 200-second- 
burning- time case. The 400 -second case is more tolerant of displacement 
errors in the lateral plane and performs comparably for the velocity 
offsets. The comparative percentages of additional characteristic 
velocity required to correct a 2° velocity offset are 18.3 percent for 
the 200-second case, and 21.3 percent for the 400 -second case, although 
the additional magnitude remains about the same. For the worst lateral 
displacement situation, ±20,000 feet about the nominal for the 2° veloc- 
ity offset, the percentages of characteristic velocity which must be 
added are 6. 5 percent for the 200- second case and 3*6 percent for the 
^-00- second case. 


Simplified Guidance Equations 

A simplified guidance technique was also tested to determine the 
importance of the corrective terms retained in the normal guidance 
relations. The only changes involved the equations for computing T x /m, 

the Z and Z commands, and for generating the measured Z and Z 
from the radar data. Equations (l 5 )> ( ) , and ( 17 ) were simplified to 



Z-com = -P (A3 + A^p^ 

and ^ttieas and ^meas from appendix E were simplified to 

^meas = P a 
^meas = P a + P^ 

Only the nominal case and the extreme errors were tested by use of 
this system. Table II is a comparison of performance for the normal 
and simplified cases. For most of the coplanar cases, a small increase 
in fuel consumption was required by the simplified system for large 
errors. For the off set- velocity cases, the simplified system appeared 
to be slightly more economical for large errors. 
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Comparative evaluation of these two systems leads to the conclu- 
sion that in most cases the simplified system would be indicated. . For 
large vehicles and particular error probabilities, the increased instru 
mentation weight of the normal system might be indicated. 


Overall Mass-Ratio Determination 

The variation of overall mass ratio (fuel consumption) with burning 
time of the terminal stage can be determined only by looking at the . 
entire energy requirement for launch and rendezvous. In order to sim- 
plify this evaluation, an impulsive launch was considered. The constraint 
applied to the comparison was that the transfer angle from launch to the 
initiation of terminal- stage burning should be the same for each burning 
time . 


Figure 12(a) is a plot of the variation in characteristic-velocity 
requirements at launch and at rendezvous as a function of the angle of 
transfer from launch to the commencement of terminal guidance for the 
coplanar case with terminal- stage burning times ranging from 100 to 
MX) seconds and for the two-impulse case. Empirical relations were devel- 
oped from the previously described exact calculations and used with 
Keplerian orbital relationships to compute the required characteristic 
velocities. Details of this work appear in appendix A. 

These characteristic-velocity requirements are not significant 
without esti ma tes of the efficiency of the launch and terminal stages. 
Accordingly, a comparison of the overall mass ratios was made subject 
to the assumption that the launch velocity is gained at a specific 
impulse of 270 seconds and the terminal velocity is gained at a specific 
impulse of 311 seconds. The relation used to calculate mass ratio was 

/ V o AV \ 

VA=e \ gI sp,l gI sp,2/ 

The variation of mass ratio with transfer angle is shown in figure 12(b). 
As expected, the two-impulse transfer is the most efficient. However, 
this is an idealized case and practical considerations will require a 
burning t ime of some length, particularly for the high te rminal- stage . 
velocity gains dictated by the optimum mass-ratio transfer angles indi- 
cated in the figure for the impulsive case and the shorter burning times. 
The percentage increase from the optimum two-impulse transfer to the 
optimum 4-00- second terminal- stage -burning transfer is the order of 
1 percent. In the real case the best burning time would be the minimum 
capable of correcting the errors expected in launch guidance. When 
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compared, with the shortest burning time practical, the penalty asso- 
ciated with adding time to correct errors will generally be much less 
than 1 percent. 

The variation of mass ratio with lateral offset was determined at 
the 82 . 5 ° transfer angle studied for the impulsive, the 200-second- 
burning-time, and the 400-second-burning-time cases. The 200- and 
400-second cases were run in the rectangular axes discussed previously 
since they were controlled. The characteristic velocities derived from 
the computations for these two cases were amplified by the ratio of 
characteristic velocity required for the exact nominal case to that 
required for the nominal case in rectangular axes to make their magni- 
tudes significantly comparable to the impulsive case. The linearity 
of the plot of the ratio of characteristic velocity required against 
initial-circumferential-velocity difference in figure 5 indicates that 
this is a valid weighting technique. 

The mass ratios resulting from this calculation are plotted in 
figure 13. The impulsive case is again the most efficient, and it is 
seen that the mass-ratio penalty for longer burning times increases 
with the velocity offset angles. 


CONCLUDING REMARKS 


The intent of this paper is to describe an automatic terminal 
guidance concept for satellite rendezvous and to illustrate some of the 
steps necessary to implement the system. The guidance system described 
has a nominal mode of operation, using constant thrust in the local hori- 
zontal plane to close the velocity difference in the terminal stage of 
direct-ascent rendezvous with an orbiting space station. Initial con- 
dition errors for the terminal stage will generally occur as a result of 
launch guidance inaccuracies. The system nullifies these errors by 
varying the magnitude and direction of the thrust according to commands 
calculated from onboard measurements of the relative position and veloc- 
ity of the vehicles. A control scheme for tracking these commands was 
devised and its dynamics analyzed. 

Sample calculations were made for a particular case, rendezvous 
with a satellite in 400- nautical- mile circular orbit for an 82 . 5 ° trans- 
fer angle. The range of initial errors with respect to a nominal aim 
point which the system could correct was determined. The limiting 
factor was system stability. For a 200-second terminal- stage burning 
time, initial circumferential errors of ± 25,000 feet in displacement 
and 300 to -400 ft /sec in velocity, and radial errors of 7,000 to 
-9,000 feet in displacement and 180 to -200 ft /sec in velocity could 
be corrected. Lateral velocities were tested only for offset angles 
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up to 2°, and it was found that displacements of ±20,000 feet from the 
nominal aim point could be corrected for all cases. The 400-second 
case demonstrated capability to correct larger errors. Provision of 
about 200 ft/sec additional terminal- stage characteristic-velocity 
capability (15 percent additional) was sufficient to handle all the 
in-plane errors if the initial circumferential velocity is not low by 
more than 150 ft /sec, and out-of-plane errors up to 1 velocity offset. 

A more sophisticated control system should be capable of correcting 
larger errors and should demonstrate greater economy. A simplified 
set of guidance commands was tested and worked successfully with vir- 
tually no additional fuel consumption. 

General analysis of the coplanar launch and rendezvous energy 
requirements for this case led to the conclusion that the shortest 
burning time consistent with the size of errors expected from the launch 
guidance should be used for the best mass ratio, and demonstrated that 
the penalty in mass ratio for using low thrust and burning times longer 
than the minimum possible was less than 1 percent. 


Langley Research Center, 

National Aeronautics and Space Adminstration, 
Langley Field, Va. , May 19> 1961- 
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APPENDIX A 

EXACT CALCULATIONS FOR FINITE BURNING TIMES AND 
IMPULSIVE TRANSFER VELOCITY - REQUIREMENTS 


In order to determine the degree of approximation involved in using 
the simplified equations of motion, trajectories were numerically cal- 
culated on the IBM type 704 electronic data processing machine for the 
terminal stage by using the exact dynamics. The exact equations were 
written in polar coordinates. These cases were for constant thrust in 
the circumferential direction, and were run in negative time from the 
rendezvous condition for burning times of 100, 200, 300, and 4-00 seconds. 
Negative thrust was used to drive the ferry to the initial condition 
from which normal rendezvous could be accomplished. The equations pro- 
gramed were 

r = R(n ) 2 - 


s - 0 - 2 ™) 


mm t 

m = “O " ~ t 

X = R S Q 
X = R s n 

Z = R s - R 

• • 

Z = -R 

This program was used for a range of thrusts to determine the variation 
of terminal- stage characteristic velocity with initial circumferential 
velocity increment to be gained for rendezvous with a station in a 
400- nautical-mile circular orbit. Characteristic velocities based on 
the mass change during these runs were calculated by means of equation (13) 
Ratios of these characteristic velocities to approximate characteristic 
velocities were calculated and are plotted against the initial circum- 
ferential velocity increment in figure 5* 

In order to determine the loss of efficiency with increase in burning 
time, the initial launch velocity impulse required to establish the 
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desired terminal- stage initial conditions was calculated where vacuum 
trajectories were assumed. Empirical relations were developed for the 
terminal- stage initial conditions from the program described previously. 
The constraint placed on the launching trajectories for comparison of 
the various burning times was that the angle of transfer from launch to 
ignition of the terminal stage should be constant. This constraint was 
felt to be reasonable, inasmuch as it represented a condition where launch 
trajectory errors would be the same for each case at initiation. Com- 
putations were carried out for assumed circumferential velocity increments 
at te rmin al- stage ignition and the results cross-plotted to yield the 
variation of required launch velocity impulse and of final circumferen- 
tial velocity increment at ignition with transfer angle. These results 
are plotted in figure 12(a). The relations used in obtaining these 
quantities, with V s - given, are (from definitions given in the 

following sketch) : 
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K± - Rg 


V S (V S - Vt,i)( At ) 2 

2.95H S 



Li 

L 

3 

d 

:> 



The expressions for V R ^ ± and R* are empirical, and the remainder are 
derived from elliptical-orbit equations. 
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APPENDIX B 


APPLICATION OF WRIGLEY ’S PROPORTIONAL NAVIGATION 
TO STEERING IN THE HORIZONTAL PLANE 


In reference 5, Wrigley proposes a proportional steering scheme 
whereby the thrust is applied in such a manner that the angular rate of 
the line of sight is driven to zero. The specification for the thrust 
acceleration vector is 


fc - S cs o> Zs x v c , s + i v v c>s + (o^ s terms) (Bl) 


where 

f c thrust per unit mass of commuter 

S cs proportionality constant 

a3 ls rotation of line of sight in space-station coordinates 

v „ commuter velocity relative to station 

C f s 

2 V unit vector along V c> s 

5j. s rotation of space station in earth- centered inertial 

’ coordinates 

The aS T „ terms will be neglected since they can be shown to be small 

in comparison to the other terms. If i, j y and k are defined as unit 
vectors along the projection of the line of sight in the horizontal plane, 
the normal in this plane, and the Z- station axis, then the vectors in 
equation (Bl) can he written in terms of these coordinates. Refer to the 
following sketch: 
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f c, s 


L v v 

v c, s 


= i[(v s - x) cos p - Y sin p] - J [(v s - x)sin P + Y cos pj (B4) 
[(V s - x)cos P - Y sin p] _ f(V B - x)sin P + Y cos p] 

' [(Vs - f * ^ 


r / . N 2 -o'] 1 / 2 

[(Vs - x ) + H 


YY - (V s - X)X 

' c ’ s 


(B5) 

(b6) 


and the Z-velocity and Z- displacements are neglected. The accelerations 
X and Y can be written in terms of the thrust and heading angle \Jf. 

A second set of unit vectors ex>ey along the X,Y axes is defined 
for this purpose in the following sketch, and Z-velocities and 
Z- displacements are again neglected. 
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where 

( R ) s = e^X + eyY 
(r)s = e x X + e Y Y 
R = e^-X ^ SyY + kR s 
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The m 2 x term in X is negligible in view of previous assumptions. 

I s 

Substitute the identity \|; = P + t where ty is a command tilt angle 
with respect to the line of sight as indicated in the following sketch: 



and 


T ^ ^ 

X = —(cos cos (3 - sin ty sin p) 
in 


T 

X x 

m 


sin p) 

(BT) 

sin p) 

(B8) 


Substitute equations (B?) and (B8) into equation (b6). 


V Cf s 


m 


[y cos P + (V s - x) sin [ 

^jsin ty + 

Y sin p - (V s - x)cos p] 

COS ty 

[( 

• \2 
v s - x) 

¥ Y 2 

1/2 

/ \ 


(B9) 


Substitute equations (B2), (B3), (b 4), (B5), and (B9) into equation (Bl) 
and equate components of the unit vectors to yield 


Is 


• T y ~ (V s - x) cos 3 - Y sin p 
S, R P + sin ty H L 


— sin ty 
m 


(v s - x ) 2 + i 2 

[(v 6 - x)cos 0 - Y sin Pj 

if +Y 2 


[(v s - x) 


>] 


sin p + Y cos P 


(BIO) 


(v s - x 
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T s - *)» 

(vT 7 

2t sln , [( Y S - x) 5 ln P ♦ i COS pj~ 

" (v s - x ) 2 + i 2 

The longitudinal guidance system will command the value for T x /m. 

Tilt the thrust vector away from the line of sight by the angle \j? to 
generate the required and assume 



f n = 


S csP + cos * “ 


n 3 + Y cos 3 
X) 2 + Y 2 


- [(v s - x) 


XJcos p - Y sin 3 


sin ^ ^ cos « 1 


then 


_ J-Y (%, 

m 


and 


*>+ 1 

♦ =-^- 

m *- 


, ^ . T x (v s - x)sin 3 + Y cos 3 

5 csP + — , . \2 • P 


( V S - it) + i 2 


(v s - x) cos 3 - Y sin 3 


[(v s - x)sin 3 + Y cos 3 


(v s - x ) 2 + Y 2 


(B12) 


Now introduce the following expressions derived from equations (10) 
and (C6): 


x 

m 


(y s - *f + * 2 

2 [( X - Vs ^ 2 + Y 2 


1/2 


(B13) 


3[(x - V s t) 2 + Y 2 ] 1 ^ = (v s - x)sin 3 + Y cos 3 
Substituting these expressions into equation (B12) yields 


(B14) 


* = 


{^cs + 


cos 3 - Y sin 


m 


l)[( V B - *) 

x + (X - V s t) 2 + Y g 

(v s - kf + i 2 


0] 


(B15) 
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L 

1 

5 

2 

2 


2 

For the conditions considered in this report, the & term in the 
denominator is negligible. Moreover the guidance system drives (3 
toward zero. Neglect this term and define 



+ 


1 

2 


to yield the final vjr expression, which is 


* = 




Y sin 


»] 


(Bl6) 


Now inspect the thrust required along the line of sight^ Introducing 
equations (Bl4) and (Bl6) into equation (BIO) and assuming \|r to be small, 
as before, yields 



Again neglect 3 ; the following result is apparent: 



and the required thrust is directed toward the space station. Accord- 
ingly, no lateral correction was introduced and the longitudinal system 
was allowed to control It/m. 
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APPENDIX C 

INTRODUCTION OF MEASURED VARIABLES INTO 
GUIDANCE COMMAND EQUATIONS 


Express the variables measured by the rendezvous guidance system 
in terms of the rectangular- coordinate system geometry. (See fig. 5-) 


P = 



(Cl) 


. _ (X - V s t)(X - V s ) + YY + ZZ 
P P 


(C2) 


I 

3 


a = sin-l(|j 

(C3) 

(p 2 - z 2 ) l/2 

(C4) 

" ■ tan_1 (x -\.t) 

(C5) 

y(x - v s t) - y(x - v s ) 
(x - Vot) 2 + Y 2 

(c6) 


For simplicity, consider that each measurement is made with the station 
at the origin of a new set of coordinates and t = 0, and use the meas- 
ured variables with a subscript 0 to denote this. 

It is necessary to restate equations ( 10) to (12) in terms of the 
variables sensed by the rendezvous guidance system. Express V s - £ 0 

and X Q in terms of the measured variables p, £>, and a. The quanti- 
ties Y and IT must be neglected, since it is not possible to isolate 
these quantities in the proposed measurement scheme. 



3K 


33 


Vc; - Xq = pa sin a - p cos a 


(C7) 


X n = p cos a 


(08) 


Represent the sine and cosine functions by their respective senes. 
Since maximum values of elevation angle are expected to be less than 
0.3 radian, truncate these series after the cubic power term 


sin a 


a - 




(C9) 


~ i a 

cos a « 1 - — 


(CIO) 


Note that the sine and cosine terms can be retained and supplied to the 
guidance computer by a resolver. It is felt that supplying the eleva- 
UofTgle d^rectl/vould result In both a veight saving and an increased 
reliability. Two approximations can be used to help weight terms 
arriving at final expressions. These approximations are 



(Cll) 


(C12) 


Now equations (C7) to (C12) can be substituted into equs.tl.ons ( 10) , 
and (12). Retaining terms consistent with previous approximations 


( 11 ), 

yields 


^com 




(C13) 


Z>com = ~P 


( 1 ■ t) 


a 4 p 


(Cl4) 



where 



5 ^ 


and 


A 2 


5H S 


A 


3 





(C15) 


Use similar means to introduce measured variables into the lateral- 
control command equation, which is 


t 


com 


= *4" [( V s - Xo)cos P - Y 0 sin p] 


J2F 

m 


(Cl6) 


Note that sin 0 O , cos 0 O , and p can be expressed 


sin p 0 = 


cos p o = 


p cos a 


p cos a 


• (*° ~ V s) X o + Y o*o + Z o*o 


(ciT) 

(Cl8) 

(C19) 


Substituting equations (C17) to (C19) into equation (Cl6) gives 


r com 


KaP 

= -2-(p 

T. 


cos a - p& sin a) 


_2L 

m 


(C20) 


Now the same means of weighting terms used before can be employed with 
the additional consideration that the lateral and longitudinal systems 


t-3 H lA CVJ CVJ 
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should be divorced as much as possible, and the rather simple control 
relationship results 

K S^ 

*com = t & 
m 


(C21) 
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APPENDIX D 

STABILITY ANALYSIS OF LONGITUDINAL AND LATERAL CONTROLS 


Longitudinal Controller 

In this appendix, the notation and techniques will be drawn 
principally from reference 10. By referring to figure 6, write transfer 
functions for the various blocks in the diagram. If the vehicle dynamics 
are linearized 


T z ^B^y 


4 


M y ^x^8®y 

Wy 

I V 


T T 

Z = — cos 0 sin 0 

m m 



m m 


Now taking the Laplace transforms of these variables and writing them 
in transfer- function form yields 

q_ ^x K 6 1 

By"' Iy X 
J_ = _ _1_ 

6y J y X 2 

Z _ ^x ^x^B 1 \ 1 

6y " V m ® Iy X 2 /x 2 


^y \ d x T x 



^ H lT\ OJ OJ 
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Consider the attitude control loop shown in the following diagram: 



The open- loop performance function for the inner loop is 

We - 4 ¥)* s < 

S p S q d x K 6 1 

I y A 


Let 


T 


_i2_ 


then 



Closing this loop gives 



1 


1 _ 

tA 


+ 


T A 


= 1 
tA + 1 

and writing the overall loop transfer function gives 


* 

A SqdxKB t A + 1 
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The open-loop performance function for the outer loop is 



a ( VeW s 

Sqd x K5 tA + 1 ^ Iy / ^2 


Sq 1 
Sq A(tA + l) 


Closing this loop gives 



S 8 1 
S q A(tA + 1) 

1 + £e 1 

Sq a(tA + l) 


— ^ A 2 

s e 



i 


i 



where 




2£ _ fq 
® s 9 


and writing the total attitude loop transfer function gives 



I 


z 


S0dxK & 



It is apparent that this simplified inner loop is of second order and 
has dynamics which may be arbitrarily selected by the designer. 
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Now consider the total control loop illustrated in the following 
diagram: 



The Z com can be combined with Z com , with the reasoning that the non- 

linearities in this quantity are gravitational in origin and may be 
neglected if gravity is neglected. In fact, these commands serve to 
make the control system detect a nearly gravity-free environment and, 
indeed, to make the linearized analysis a good representation of the 
system's real performance. The resulting block diagram is 



The open- loop performance function for the inner loop of the 
preceding sketch is 




Closing this loop will yield a third-order characteristic equation. 
The result will be of the form 
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where t^, and can be determined by root- locus techniques. 

Again this function can be carried into the outer loop, and the 
open-loop performance function determined 



For the guidance-problem studies in this paper, the attitude loop 
was arbitrarily set up for (; = 0.7 and o> = 1.4l4 radians. This con- 
dition yielded the following system constants where an arbitrary value 
of 300 was selected for SpK§ and 



Note that S p is redundant. 

The root locus for the inner Z-loop appears in figure 14 with the 
closed- loop roots selected for the 200- and 400-second burning times. 

These roots were carried to the outer loop plot, and the two corresponding 
loci also appear in figure l4. The final Z-loop gains corresponding to 
the open- loop sensitivities from the root loci were as follows: 


400- second 
case 

0. 18 
0.008 

Note that the response for these gains is rather slow. These gains 
were purposely selected to allow the system to operate on very large 
errors before exceeding the tilt-angle limit, and since a relatively 
long time is available for correcting errors. 


KZ 


200-second 

case 

0.122 

0. 0091 



H LT\ OJ OJ 
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Lateral Controller 

Refer to figure 6 and note that^in this section the object is to 
establish the commanded tilt angle \j?. Linearizing the lateral dynamic 
equations yields 

Ty = K 6 5 z 
M z = d x K S & z 


M„ 


K*d. 


r = 


5 a x 


Taking the Laplace transform yields 


£ K s d x i 

^z Iz Tv 2 


Consider the inner loop shown in the following sketch 



and write the open-loop performance function 



_ _ K 6 d x 1 , a 
~ S L ~ f — ~2 

-L?. A 


SlWx 1 

Iz A 


where 


t 2 = 


I 


z 
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Then 


If this loop is closed 



1 

t 2 A 


[pf] 


cl ToA + 1 


and the desired transfer function can he written 


i = _i_ i 

^ Sj»A t^A ■+■ 1 

Now follow the same procedure for the outer loop to obtain the fol- 
lowing equations: 


Nor 

t = 

*com 

where 



2^2 fr 

^ = 


The dynamics of this system are again at the disposal of the designer. 
The following gains, corresponding to £ = 0.7 and = 1.428, were 

selected with S-^ being redundant: 

s l k & = 300 

s* = 5 


_■ ft 1 

S r A(t 2 A + 1) 


■A.\ a + 54 + i 

“2 


I 05, 


S,S T K e d. 


S r = 2.5 


hJ H 1A(AJ OJ 
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APPENDIX E 

EQUATIONS OF MOTION FOR NUMERICAL SOLUTION 


The equations of motion for the rendezvous vehicle were written 
in Eulerian axes centered at the vehicle center of gravity. (See fig. 1. ) 
Since the vehicle was assumed to be roll stabilized to the horizontal 
plane sensed by the horizon seeker, the third Euler angle 0 rotation 
was suppressed. This has the effect of constraining the Y-axis to the 
horizontal plane and requires the following substitution to account for 
rotation about the X-axis as a result of this constraint: 

0 = p - r tan 9=0 


p = r tan 0 


The motion of the vehicle was referred to a rectangular X, Y, Z axis 
system with origin at the space station's location at time zero. Addi- 
tional auxiliary relations were required to calculate the measured 
variables and to represent the control system; these relations are 
listed in appropriate groups. 


Dynamics: 


T 

u = rv - qw + — - G sin 0 
m 

T 

v = -r(u + w tan 0) + -£ 


w 


rv tan 0 + qu + — + G cos 0 
m 


4 


r 



6 = q 


if 

G 


r 

cos 0 


(Rs - Z)2 


X^Y 2 
R s - 


Z 



kk 


Varying mass parameters: 


m = - -**■ 
c 

d x = ^,0 + d x,l t 
J x = + I x > l t 


I y " h,0 + I yA t + I y> 2 1 " + -’y, 3 1 


I z - I z ,0 + I z,l t + I z ,2 t + I z,3 i 
Geometrical relations: 

X = u cos \|r cos 0 - v sin \|r + w cos i|r 
Y = u sin >|r cos 0 + v cos i|r + w sin 
Z = -u sin 0 + w cos 0 
Measured variables: 

p = [(X - V s t) 2 + Y 2 + Z 2 J 


. 21 1 / 2 


= (X - V s t) (X - V s ) + YY + 


zz 


a = sin 


a - 


1 !) 

1 - 4 ) 

(p 2 - Z2) l/a 



. r(x - v s t) - y(x - v s ) 

" 2 o 

(X - V s tr + ?r 

? = t - 3 


3 

3 

sin 6 
sin 9 


ni h incvj oj 
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Control forces and moments: 





m 



Mz = K 6 d x 
Iz d Z 


5 


z 


Guidance commands : 



Feedback variables: 

^meas 

Zmeas 



■ 4 - ir) 

• 4 - t) * 4 • 


Fuel consumption relations: 
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* - 

where | j denotes absolute value. 

The physical constants employed were derived by scaling the Mercury 
capsule up to a 300-slug vehicle and elongating it to about 18 feet. 

The control constants appear in appendix D, and the remaining physical 
constants are tabulated as follows: 


R s , ft ... . 
V s , ft/sec 2 . 

7 , .ft3/sec 2 
IDq, slugs . . 

l x , 0> slug- ft 2 

ly, 0> SlUg- ft 2 

l z, 0 > slug- ft 2 

d x, 0 > • • • 

c, ft/sec . . 


• • • 23 , 333,580 

. . . 24,563 .66 

1.4078741 x 10 16 

300 

1,800 

7, 500 

...... 7,500 

10 

10, 000 


200 - second 
case 


400 - second 
case 


^x, 1* 

ft/sec 


-0.0025 

J x,l> 

slug-ft 2 /sec 

. -1.225 

-0.6125 

T y, i' 

slug-ft^/sec 

- 7.5 

-3.75 


slug-ft 2 /sec 

- 7.5 

- 3-75 

x y, 2 > 

slug-ft 2 /sec 2 

- 0.5 X 10“ 2 

-O.125 x 10‘ 2 

Iz, 2 > 

slug-ft 2 /sec 2 

-0.5 X 10 “ 2 

-0.125 x 10~ 2 

T y, 3> 

slug-ft^/sec^ 


-0.625 x 10~6 


slug-ft 2 /sec^ 

-0.5 x 10" 5 

-0.625 x 10“ 6 
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TABLE I.- NOMINAL INITIAL CONDITIONS 


(a) 200 - second burning time 


L 


o, deg 

0 

1 

2 



135,425 

129,525 

115 , 84 i 

y 0 , ft 

0 

39, 537 

70, 149 

z 0 , ^ 

18,960 

18, 097 

16, 094 

Uq, ft/sec 

22,239 

22, 104 

19 , 446 

v Q , ft/sec 

0 

- 7,171 

- 12,492 

w Q , ft/sec 

-280 

-280 

-280 

radians 

0 

0.2963 

0. 5445 


(b) 400 - second burning time 


0, deg 

0 

1 

2 

ft 

239,619 

226, 424 

197,563 

Y 0 , ft 

0 

78, 420 

135,596 

z 0 > ft 

67,214 

63,000 

54, 553 


23,389 

21, 964 

18, 810 

v Q , ft/sec 

0 

- 8, 039 

-13,900 

w Q , ft/sec 

-496 

-495 

-496 

ty Q , radians 

0 

0.3334 

0. 6015 


ro ro vji h 
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TABLE II.- COMPARISON OF REQUIRED CHARACTERISTIC VELOCITY 
FOR NORMAL AND SIMPLIFIED GUIDANCE RELATIONS 


Characteristic velocity, AV 


Error 


For normal 
guidance relations 

a = 0° 


Nominal 

1, 325.00 

AX = 25,000 

1, 499.48 

AX = -25,000 

1, 471.23 

AZ = 7,000 

1, 498.31 

AZ = -7,000 

1,484.41 

AX = -400 

1,842.02 

AX = 300 

l, 156.99 

a£ = 180 

1, 508.88 

AZ = -200 

1, 504.79 

a = * 

! 

Nominal 

1,389. ^5 

AY = -20,000 

l, 445.68 

AY = 20,000 

l, 454.63 

a = i 

Nominal 

1,566.82 

AY = -20,000 

1,647.57 

AY = 20,000 

l, 667.33 


For simplified 
guidance relations 


1, 325.55 
1 , 505.23 
1,475-77 
1, 502.82 
1,491.89 
1,844. 51 

1.163.67 

1. 507.67 
1, 505.01 


1, 389.72 
1,438.60 

1, 452.21 


1,567.75 
1,658. 38 
1, 656. 05 
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Figure k.- Comparison of rendezvous displacement and relative velocity profiles from exact cal- 
culations and assumed rectangular representation for matched initial circumferential veloc- 
ities with 200- and ^00-second terminal- stage burning times. 





Figure k m - Concluded. 
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(a) Relative geometry and thrust variation. 

Figure 7 .- Coplanar rendezvous for nominal and extreme initial errors 
in displacement. Nominal terminal- stage burning time of 200 seconds 
for rendezvous with station in 400-nautical-mile circular orbit. 
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(b) Time history of final 10 seconds before rendezvous. 


Figure 7* - Concluded, 
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(b) Time history of final 10 seconds before rendezvous 
Figure 8 . - Concluded. 
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Figure 10.- Time history of final 10 seconds before rendezvous for ini- 
tial lateral velocity offset of 2° at nominal initial lateral offset, 
extreme initial displacement errors, and Kp = 3«0* Nominal 

terminal.- stage burning time of 200 seconds for rendezvous with sta- 
tion in 400- nautical-mile circular orbit. 
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200-second burning time 0 400-second burning time 



(a) Initial in-plane errors. 

Figure 11.- Variation of required characteristic velocity with initial in-plane errors and cross 
plane displacements for burning times of 200 and 400 seconds. 
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(a) Required characteristic velocities. 

Figure 12.- Characteristic velocity and variation of mass ratio for 
two-impulse coplanar transfer and finite terminal- stage burning 
times of 100, 200, 500 , and 400 seconds. 
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Figure 12.- Concluded. 
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Figure 13 .- Mass ratio variation with lateral velocity offset angle for two-impulse case 

terminal- stage burning times of 200 and kOO seconds. 


68 








u ~ u m 

T3 S3 - 


73 3 B § 
■2 O <i> n 

S - » S 

bb cs co b 


^ S « 
£ o qj 13 


•*2 * w m 

:H |S«fra 


Q <s 
Eag 
< § < 
«3 d Z 
zz< 


W »Z 

H|B 

ig«£ 

3 W 5*2 

< D c < 

>o sS 

; > o 

< N . 

!q 2 - 

!g 

| 2 S^ 
■ os E 3 

i fa H a 


,d 4) *"' 

p 23 an <d ' 

S - b s o . 

* ® a i2 3 

0 3 .2 c 3 

3 B -8 ? 
s|g£° 

q c 2 a® , 

1 S Si 2d 

“ a> * «« d ■ 
o S3 *0 v - 
M a u o c : 
S 3 S -2 
a> d 3 13 ; 

l*% vi w 

»13 |d 
o Si 6a5 

^3 d d co ^ 

c A 6 p ^ 

O y 5 CO . 

u S 3 ■ 

— u 2 C O 

^ .2 3 M-* w 

5 T3 £ T3 N ■ 

3 fi 5^- 

b£ .5 N TJ 

,5 C! 3 J^T* <JJ 


CO *0 
>,oj c a 

£ 3 p s 

| -2 tS 3 

> c 'O a> 

73 5 * ® 

II 8 8 


« '3 " 

JC CO u 3 


•S CD ° ^ CD 

V, CJ ^ c3 U 

a> c o c 

® 4 is 3 ©“ 

co .b e *o co 


< a 8 s s -s 


s^iS 

u w S CO 

?sl* 

£« « u * 
3 S 4> § 
3 0 > o3 'S 


jw<§-« 

fc^g 

J .2 

Jh , 

£S« 

; 5 
H §>g 

il<g 

sr| 

* D c 5 
)? &£ 

!gs« 

> £ <u 2 . 

>s § . 
: ^ 2 « 
g D u 
;0 ait 
ifcha 


T3 CD 

73 boi a 

3 5 g a 

g s gs*a i 

•2 g a 3 a' 

a 3 .2 c g 

n > s3 

a | g g § ; 

im:- 

* -S c ■ 

o JS U . 

^ 3 T3 O C ; 

a ^ g ^ . 

£3 td , 
W QO 0) -5 

* ^ T3 ti Cfl 

*f 5 |- 

2^53,5 

iJ d d yi > 

C ' H -3 

O y _ 5 OT ' 

iSlag: 

bO -3 n ^ 

.S <d O >» a) 

t( u h d t 

oj c o c 


3 g CD 

5 o 0) d 

D -2 O CO 

1 1 a- 

■d 3 m - 


^ -2 O CO 

«■« a . 


S8 « 
z | n I 
.3 o S 5 
•2 

£ <n >» S 

OCN3 CO b 


dJ 3 £ 5 

.3 o S 5 

•3 - to o 

£s s-p 


a w® g 

J .2 

Jh , 

3«s< 

: w a s 
H|g 

5 £ < a 

b n 


T3 cu 

73 MS a 

3 s g a 

3 2 CO CD ' 

2 d V > 13 * 

a; ~ b <«-i CJ i 

® ® a 3 xj 

u 3 .2 c 3 

d N " M 7 
E^oSg 
q c a® , 

d ^ g ^ d 

U p i rt S . 

O .tS T 3 Vi ■ 

^Ztj.oc 

b 0 ) fl &M O 

§ d ™ 13 ; 

ti » 0 ) ■ 2 

* t) tn 03 

»|3 a- 

oS§3 a5 
is rt d co *> 

c ' gd 

S^U5 5 3- 

H .§ 3 -S § 

S -3 -a N • 

3 c £^' 

■ S CD O 5* CD 

in cj — < d v 

<d c o d 

CO -b d TJ 50 

. 3 o c d . 

< & y d y 


^ O . 

*2 w S 

5 o w 

d _< 

co T3 
>> v d 2 

PH 

-r! CO Vi d 
g - <D X 
> g -a a) 

73 S W <D 

3 " e H 

ti co o d 

• S Cfi b 


® c g 3 
b 2 v, -g 

iS £ « S 

3 Qf q T 

y y -o b 


Ih 

eo ^ y 
« c r 

ai o cl l 
i u <c | 
0^5 
< S 

Z Z P 

^ gSi 
|-2<, 
< d Z < 

zz^l 


"O ni 

73 bDii a 

I .S g E 

g s g | -a 1 

" « a 3 3 ' 

s i " 8 I 

g g c 2 • : 

I I -2 s.8 . 


o bl T3 V . 
41-1 a "d o d 1 
C 4) C .2 

a) d ™ 13 ; 

tb co (D .5 

2 _ T3 V CO 

» |3 
o g aSs 

d d co > 

it ss 
" £ » = S 
S^r^s: 


XJ *-* v- 
50 -b C -3 
^ 3 0 d rs - 





